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We investigate the low-lying excitation of the antiferromagnetic zigzag spin ladder with Mobius 
boundary condition. Using the Lanczos and Householder diagonalization methods, we calculate the 
excitation spectrum of the zigzag ladder in the momentum space. We then show that the topological 
defect generated by the Mobius boundary provides the clear evidence of the single spinon excitation. 
On the basis of the obtained single-spinon dispersion curve, we analyze the two-spinon scattering 
states under the usual periodic boundary condition. We further discuss the interaction effect between 
the spinous, and the connection to the cusp singularity in the magnetization process. 



Antiferromagnetic(AF) S — 1/2 Heisenberg spin lad- 
der with the zigzag structure has been an important issue 
in theoretical studies of quasi one-dimensional(lD) quan- 
tum spin systems, since it is a minimal model including 
both of the quantum fluctuation and the frustrated in- 
teraction. j^J^ Particularly it is an interesting problem 
to clarify how the effect of the frustration appears in 
the structure of the low- lying excitations, which is 
closely related to various characteristic behaviors of ob- 
servable physical quantities, such as the cusp singularity 
in the magnetization curve. [^] 

In analyzing the low-energy excitation of the zigzag 
antiferromagnet, an essential object is the spinon excita- 
tion, which may be described by a "local-defect particle" 
in the ground state. However, the nature of the single 
spinon excitation can not been investigated directly for 
the system with the periodic boundary condition(PBC), 
which is assumed in the most of theoretical studies; The 
total-S* of the system under the PBC takes integer values 
and hence one can not extract the single-spinon excita- 
tion of 5 = 1/2 straightforwardly. 

In order to discuss the topological-defect excitation, 
in this letter, we deal with the AF zigzag spin ladder 
with the Mobius boundary condition(MBC), which cre- 
ates a domain-wall type defect in the system without los- 
ing the translational invariance.(See Fig.|l|) In contrast to 
the usual PBC, the zigzag ladder with the Mobios bound- 
ary contains the odd number of spins due to the network 
structure of the zigzag lattice, where the total-S* of the 
system can be labeled with a half integer. |^ Then an es- 
sential difference can be expected in the low-energy spec- 
trum for the MBC, where the defect particle is able to 
move with a finite momentum. We find that this defect 
particle, i.e. the Mobius twist in the zigzag spin ladder, 
is identified as the spinon excitation of 5' = 1/2. 

In the following, we discuss the spinon excitation gen- 
erated by the MBC, using the Lanczos and Householder 
diagonalization technique in the momentum space. On 



the basis of the calculated spinon dispersion for the 
MBC, we further consider the low-energy excitation of 
the zigzag spin ladder with the PBC, where the two 
spinon scattering state plays an important role. In addi- 
tion, we show that the shape change of the spinon disper- 
sion curve, which is connected with the cusp singularity 
in the magnetization process, can be illustrated in the 
exact diagonalization spectrum obtained for the MBC. 




FIG. 1. The antiferromagnetic zigzag spin ladder with the 
Mobius boundary. The inter-ladder couplings, which are nor- 
malized to be unity, are denoted by thin line and the in- 
tra-ladder couplings(a) are indicated by the thick line. 



The Hamiltonian of the zigzag spin ladder is given by 

(1) 



Af-l 

7i = ^ ^Sf Si+i + aSi ■ Si+2 



where Si is the 5* = 1/2 spin operator at i-th site and a 
is the ratio of the inter- and intra-ladder couplings. The 
number of sites is denoted as TV, which takes odd integer 
for the MBC and even number for the PBC. Here we note 
that, for the infinite-length system, the ground state of 
Eq. (|^) for a < ac = 0.2411 is in the spin liquid phase 
and that for a > '\s in the dimer gapped phase. 

An important feature of the Mobius boundary is that 
the translational invariance of the system is maintained, 
although the domain wall is inserted in the system. This 
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implies that the topological defect can move around in 
the "Mobius ribbon" smoothly. (In the Fig. |l], the 
twisted part can move in the whole region of the "Mobius 
ribbon".) Thus we can make the Hamiltonian (|]) block- 
diagonal with respect to the momentum quantum num- 
ber ki = 2tiI/N, where / = 0, 1, ■ • • , — 1. We perform 
the Lanczos diagonalization up to 28 sites and the House- 
holder diagonalization up to 18 sites in the momentum 
space. 

Before proceeding to a discussion of the computed re- 
sults, here, it should be noted that, in the spin liquid 
phase, there are gapless excitations near the ground state, 
which requires us to analyze the calculated data precisely. 
Thus, in this paper, we consider the dimer-gapped phase 
ttc < a < 1 mainly, although the argument of the MBC is 
also available for the spin liquid phase. The single-spinon 
excitation in the dimer phase has been discussed only 
near the Majumdar- Ghosh point(a = 1/2), using the 
variational approach combined with the matrix-product 
ground-state. In the context of the spin-Peierls sys- 
tem, the dynamical properties of the zigzag spin ladder 
is also investigated with a variety of the numerical tech- 
niques. [^,0 However, the topological boundary condi- 
tion has not been emphasized so far. In the following, 
the striking evidence of the single spinon excitation is 
demonstrated in the exact diagonalization spectrum of 
the system with the MBC. 




FIG. 2. The one-spinon dispersion curve of the zigzag chain 
with the Mobius boundary condition, (a) Lanczos results up 
to 27 sites for a = 0.5(Majumdar-Ghosh point), 0.6 and 0.7. 
The dashed line is the spinon dispersion curve obtained with 
the variational calculation for a — 0.5. (b) The spectrum 
obtained with Householder diagonalization for a — 0.6 

In Fig. §-(a), we show the low-energy spectrum of 
the zigzag spin ladder of a = 0.5, 0.6 and 0.7, which 
are obtained with the Lanczos diagonalization up to 27 
spins in the total-S'^ = 1/2 subspace. We also show the 
Householder diagonalization results for a = 0.6 in Fig. 
H-(b), which include the higher-energy spectrum. Since 
the defect particle of the MBC can move in the ring, 
the low-energy spectrum is described with the continuous 
curve, which is adiabatically connected to the ground 
state of the MBC system. In other word, the ground- 
state energy for the MBC is given by the bottom of the 
dispersion curve of the defect particle, which is located 
at near k = Tr/2. In Fig. however, we have determined 
the origin of the energy axis with 



MBC 




(N) 



X iV, 



(2) 



where eg denotes the ground-state energy per site ob- 
tained from the ground-state energy for the PBC system 
of 28 sites with use of Cg = E^^'^{N)/N. For exam- 
ple, we have eg = -0.375(exact), -0.3802, and -0.3972 
for a = 0.5, 0.6 and 0.7 respectively, where the finite- 
size correction can be negligible. Then, we find that 
the computed spectra for the various N get on a sin- 
gle curve, which can be regarded as the one-particle dis- 
persion curve above the "virtual ground state" with the 
"virtual excitation gap" . ||ll| Moreover, in Fig. |-(b), the 
dispersion curve of the defect particle is clearly separated 
below the continuum spectrum, which implies that the 
particle is a well defined object in the zigzag spin ladder. 
Thus it will be allowed for us to call this defect particle 
in the MBC system as a "spinon excitation". In fact, 
the dispersion curve for a — 0.5 is in good agreement 
with the spinon dispersion obtained with the variational 
calculation Q , which is represented as the dashed line in 
the Fig. &(a). 

In Fig.^, we can also see that the shape of the disper- 
sion curve changes into the double well structure, which is 
expected from the cusp singularity in the magnetization 
process. Such a clear evidence of the shape change is dif- 
ficult to find for the PBC case, where the two-spinon con- 
tinuum should be dealt with numerically. (See also Fig.^) 
Precise calculations indicate that this shape change oc- 
curs at ttcusp — 0.54, which is slightly larger value than 
the Majumdar-Ghosh point a — 1/2. We note that the 
similar behavior of the dispersion curve is reported in the 
5 = 1 bilinear-biquadratic spin chain, JT^ , p^ which ex- 
hibits the magnetization cusp and the exact matrix 
product state called valence-bond-solid. [Q 
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FIG. 3. Hoeseholder diagonalization results for a = 0.6 
with the periodic boundary condition. The shaded region 
represents the two-spinon continuum. 

Let us next investigate the excitation spectrum of the 
PBC case, exploiting the spinon dispersion curve ex- 
tracted from the MBC. In Fig. we display the energy 
spectrum of the PBC for a = 0.6, which is computed 
with Householder diagonalization method up to 18 sites. 
For the case of the PBC, the dimerized ground-state is 
located at fc = and tt, and there is the excitation gap 
whose magnitude is twice of the virtual gap for the MBC. 
This magnitude of the gap is consistent with the result 
of the density matrix renormalization group. Further 
we can see that the singlet and triplet excitations con- 
form the continuum spectrum above the excitation gap, 
which can be represented as the scattering state of the 
two spinous. (In Fig.l, this will be illustrated as the two 
twists running in the zigzag ribbon independently.) 

To see this, we compare the free two-spinon continuum 
with the Householder results. The two-spinon continuum 
is given by 



£;((7) =e(fc+)+e(fc_). 



(3) 
(4) 



where e(fc) is the single-spinon dispersion curve, which 
can be determined from the result of the MBC system. 
We find that e(}i) is well fitted with the form: 



e(fc) = ^jA(k)B{k), 



(5) 



where = 1 -I- ai cos(2fc) -I- cos(4fc) and i?(fc) — 

h\ -1- &2 cos(2fc) -I- &3 cos(4fc) with the fitting parameters 
{a} and {6}. |^| The calculated two-spinon continuum 
is shown as the shaded region in Fig. ^. Clearly the 
lower bound of the two-spinon continuum can explain 
the diagonalization results, including the subtle structure 
above the excitation gap near fc = and vr. Hence we 
can conclude that the excited state in the PBC system is 
basically described with the two-spinon scattering state. 

As was seen in the above, the low-energy property 
of the PBC system is well explained with the free two- 
spinon state. For the purpose of the quantitative anal- 
ysis, however, the interaction effect between the spinous 



becomes important, especially in the finite-size system. 
In addition, we can see that the eigenvalues in the 5^ = 1 
subspaces deviate from those of S*^ = even in the low- 
energy region(see Fig. ||). This is because the frustrated 
interaction affects the excitation spectrum. Thus, let us 
next investigate the spinon-spinon interaction effect on 
the low-energy spectrum. For this purpose, we consider 
the dilute spinon gas model, which may be described with 
an effective Hamiltonian with two-body interaction: 



(6) 



fc.fc' 



where j,(as,fc) is the spinon creation(annihilation) op- 
erator and Vs^s' {k, fc') is the effective coupling. The index 
s denotes the z-component of the spinon excitation and 
the single-spinon dispersion curve e{k) has been deter- 
mined with Eq. (^. 

In order to discuss the interaction effect in the exact 
diagonalization spectrum, let us now concentrate on the 
two-spinon problem, for which the Schrodinger equation 
can be written down explicitly: 



[e{k+) + e{k^) - E{q)]Mk+,k-) = 
-^F,v(fc,fc')V's'(fcV,fc^), 



(7) 



where k± — q/2 ± k and q is the momentum of the cen- 
ter of mass. The spin symmetry can be classified in the 
spin singlet and triplet sectors. Assuming the repulsive 
(5-function-type interaction Vs^s' (k, k') = c(=const), the 
solution of is given by 



E{q)-eik+)-eik^)-' 



(8) 



We solve Eq. (^) for a finite-size system numerically and 
compare the results with the diagonalization ones. Par- 
ticularly we consider the system of the = 28, where the 
wave number is given by fc; = 2-kI/N with / = 0, 1, ■ • • , 27. 

We first discuss the spectrum obtained with the Lanc- 
zos method in the total-S'^ = subspace, where the sin- 
glet state and the S*^ = state of the triplet appear. 
Fig. ^ shows the results for the two-body problem (|^) 
with c = 0.1, where we can see the quantitative agree- 
ment in the Lanczos results of 5^ = sector. If we drop 
the interaction term(i.e. free-spinons case), the second 
lowest and the third lowest eigenvalues of the effective 
model degenerate and then a branch of the spectrum can 
not be explained. Thus we can conclude that the Lanc- 
zos spectrum of total-S'^ = is well explained by the 
effective spinon gas model with the repulsive interaction. 
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FIG. 4. Comparison of the Lanczos results with the effec- 
tive spinon gas model. 

We next proceed to the analysis of the = \ sector. 
In Fig. ^, we can see that the Lanczos results of the 
triplet sector are rather scattered, comparing with those 
of the singlet one. Such a difference between the singlet 
and triplet spectrum in the low-energy region may be a 
characteristic behavior in the strongly frustrated system, 
in contrast to the Heisenberg spin chain where the triplet 
and singlet excitations exhibit a quite good correspon- 
dence. In fact, we find that this kind of singlet-triplet 
difference emerges, as a is increased beyond a ~ 0.4. 
Accordingly, some of the spectrum of the triplet sector 
do not agree with the present analysis for a = 0.6 uti- 
lizing the simple J-function interaction; The effect of the 
frustration may appear in the triplet excitations signif- 
icantly, and the interaction term Vs^s'ik^k') can be ex- 
pected to have more complicated form than that of the 
singlet. Here, we remark that the finite-size effect in 
the triplet sector is much bigger than the singlet sector, 
which is consistent with the observation of the strongly 
affected triplet sector. In order to verify the dilute spinon 
gas model for the triplet, however, it is required to inves- 
tigate the excitation spectrum of a longer system. 

Finally we make a comment on the bound state. In 
the variational treatment of the two spinon excitation, 
Vs^s'^k^k') is not a simple repulsive interaction, where 
the bound state of the two spinon can appear near k = 
7r/2. In the diagonalization results, we can see that 

there are lower energy eigenvalues at k = 7r/2 than the 
scattering state, which is consistent with the bound state. 

To summarize, we have discussed the Mobius bound- 
ary condition(MBC) in the antiferromagnetic Heisenberg 
spin ladder with the zigzag structure. The MBC intro- 
duces the topological defect that is described with the 
domain- wall type particle of S = 1/2. Using the exact 
diagonalization method in the momentum space, we have 
found that this defect particle gives the clear evidence of 
the spinon excitation. In addition, we have shown the 
shape change of spinon dispersion curve, which is associ- 
ated with the cusp singularity of the magnetization pro- 
cess. 



We have investigated the low-energy spectrum with the 
usual periodic boundary condition(PBC), where the two 
spinon state plays an important role. We have shown 
that the two spinon scattering state explains well the 
spectrum computed with the exact diagonalization. We 
have further discussed the interaction effect between the 
spinous, using the dilute spinon gas model. The = 
sector is in good agreement with the effective model, 
while the spectrum of S'^ = 1 subspace exhibits rather 
complicated behavior. To clarify the role of the frustrated 
interaction in the triplet sector, a further analysis of the 
zigzag system will be required. 

The authors would like to thank T. Nishino for valu- 
able comments. A part of the present numerical com- 
putations was performed at Yukawa Institute Computer 
Facility. 
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